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Relativistic collisions of light on heavy ions (p+Auat y^=7.7 GeV, p+Au, d+Au ,®He+Auat 
=62.4 GeV and 200 GeV and p+Pb ,®He+Pb at y/s = 5.02 TeV) are simulated using “super- 
SONIC”, a model that includes pre-equilibrium flow, viscous hydrodynamics and a hadronic cascade 
afterburner. Even though these systems have strong gradients and only consist of at most a few 
tens of charged particles per unit rapidity, one finds evidence that a hydrodynamic description ap¬ 
plies to these systems. Based on these simulations, the presence of a triangular flow component 
in d-bAn collisions at y/s = 200 GeV is predicted to be similar in magnitude to that found in 
®He+Au collisions. Furthermore, the i> 3 (pt) ratio of ®He-|-Auto d+Auis found to be sensitive to 
the presence of pre-equilibrium flow. This would imply an experimentally accessible window into 
pre-equilibrium QGD dynamics using light-heavy ion collisions. 


2 


I. INTRODUCTION 

Experimental data on d+Au, p+Pband ^He+Au collisions seems to indicate the presence of collective flow in an 
(at least partially) equilibrated system EHI. While alternative explanation have been put forward @ , the result from 
measurements of multi-particle cumulants are strongly suggestive of a hydrodynamic origin of this collectivity Q- 

This study is motivated by the following, interrelated, questions: 

1. Does hydrodynamics apply to systems created in light-heavy ion collisions? 

2. At which system size or collision energy does hydrodynamics break down? 

3. What experimental observables would confirm or rule out a theorist’s answer to the above questions? 

In order to obtain answers for these questions, I use the new simulation package “superSONIC”, which is an event- 
by-event generalization of the SONIC model Q, including pre-equilibrium flow, relativistic viscous hydrodynamics as 
well as a hadronic cascade afterburner (details are given below). Using a Monte-Carlo Glauber model for generating 
event-by-event initial conditions, superSONIC can be used to simulate a range of different collision systems at dif¬ 
ferent collision energies, such that results are directly comparable. Comparison to data will be performed wherever 
experimental information is available. 

Quite a number of theoretical studies on light-heavy ion collisions have appeared in the recent literature, and in 
the following similarities and notable differences of these works with respect to the present study are highlighted. 

In Ref. 0, Piotr Bozek studied p-fPband d-l-Pb collisions at ^/s = 4.4 TeV and ^/s = 3.11 TeV, respectively, using 
Monte-Carlo Glauber initial conditions followed by a 3-|-ld viscous hydrodynamics evolution and a hadronic cascade 
afterburner. He successfully predicted the large flow signal in p -|- Pb collisions at = 5.02 TeV. 

In Ref. 0, Nagle et al. studied p-l-Au, d-l-Auand ^He-I-Au collisions at i/s = 200 GeV and p-l-Pb collisions at 
y/s = 5.02 TeV using Monte-Carlo Glauber initial conditions followed by a 2-I-1 viscous hydrodynamics evolution and 
a hadronic cascade afterburner (“SONIG”, the predecessor of the package used in the present study). Nagle et al. 
found results consistent with Ref. Q, and first proposed ^He-|-Au collisions as an interesting handle on QCD transport 
properties. Elliptic and triangular flow in ^He-|-Auat ^/s have since been measured experimentally [llj. 

In Refs. [H, Schenke and Venugopalan studied light-heavy ion collisions using IP-Glasma initial conditions, 
followed by a 2-|-Id viscous hydrodynamic evolution, but no hadronic cascade afterburner. Their study is similar to 
the present one in that the sudden matching from the Glasma evolution to hydrodynamics includes non-vanishing 
pre-equilibrium flow. However, given that the Glasma evolution never drives the system to equilibrium, one can expect 
the results to be strongly dependent on the Glasma-hydro switching time. In their study, Schenke and Venugopalan 
found sizable flow components V 2 , for d-l-Au ,^He-|-Aucollisions (in agreement with experiment), but very little flow 
for p-t-Pbcollisions (in disagreement with experiment). 

In Ref. [ij, Kozlov et al. studied p-l-Pb collisions at = 5.02 TeV using 3-|-Id viscous hydrodynamics, but 
without any pre-equilibrium flow or hadronic afterburner. They implemented initial conditions based on a Monte- 
Garlo Glauber model supplemented with negative binomial energy fluctuations, which turn out to be very important 
for describing rare “high-multiplicity” events. Kozlov et al. found that they could successfully describe experimental 
flow data in p-l-Pb using their model. 

Gompared to some of the studies listed above, the present study suffers from the drawback of only simulating boost- 
invariant (2-|-Id) dynamics, which clearly is not a good approximation to the longitudinal dynamics for light-heavy 
ion collisions. However, many of the conclusions in the present article will be based on ratios of flow observables in 
the transverse plane, so that one can expect conclusions to not be dominated by longitudinal artifacts. Nevertheless, 
it would be good if the present results could be checked by a full 3-|-Id calculation in the future. Also, the present 
study does not include negative binomial energy fluctuations, on the basis that these fluctuations are most relevant 
only for rare high-multiplicity events. In this article, the emphasis is on central (impact parameter b < 2 fm/c), 
but not high-multiplicity selected events. Einally, the present study is based on a geometric Monte-Garlo Glauber 
model [l^ rather than the IP-Glasma model [l^, mostly to offer a baseline result using transparent ingredients. If it 
would turn out that some experimental result could not be described by the present approach based on Monte-Garlo 
Glauber, but can be described using IP-Glasma initial conditions, this could be regarded as experimental evidence 
for Golor-Glass-Gondensate dynamics in QGD. 

On the other hand, the present simulation package for the first time combines pre-equilibrium flow, viscous hydro¬ 
dynamics and hadron cascade dynamics in an event-by-event study that is applicable to different collision geometries. 
In this sense, it is the most realistic description currently available, and the direct comparison to experimental data 
may therefore offer interesting clues about the nature of hydrodynamics and transport in strongly coupled QCD. 

The remainder of the article is structured as follows: in SecllTl a brief description of the components of superSONIC 
is given, followed by results for p-l-Au ,d-|-Au, ^He-|-Auand p-l-Pb collisions at various energies in Sec IHIl and the 
conclusions in SeclIVI 
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II. METHODOLOGY 


A. Early Stage: Initial Conditions and Pre-equlibrium 


In principle, obtaining information about real-time evolution in relativistic ion collisions is a problem that can 
be well defined in QCD. However, currently there are no known techniques to solve this problem, and thus obtain 
information about the pre-equilibrium stage in these systems. Nevertheless, a number of approximate approaches 
exits. Among these, there are weak-coupling inspired methods that rigorously apply to QCD in the asymptotically 
high energy limit (cf. [iffl), and strong-coupling methods that rigorously apply to a certain class of non-abelian 
gauge theories (but not QCD) in the limit of large number of colors and large t’Hooft coupling (cf. [26l - l3^ b 

From a modeling perspective, what is needed from any of these methods is information about when the system 
behaves approximately hydro-dynamically (equilibration time) as well as the values of the hydrodynamic degrees of 
freedom (hydrodynamic initial conditions). Other than the hydrodynamic starting time 22l. 251. results for hydrody¬ 
namic initial conditions are presently available only from the strong-coupling method [32 , 34, ^ , so this is what will 
be used in the present study. 

In more detail, from numerical relativity simulations of space-times modeling the relativistic collision of “ions” in 
J\f = A SYM it has been found that the radial fluid velocity is proportional to the gradient of the initial density 
distribution [H, HH, [s^ . In this study, this finding is promoted to the fluid velocity in all directions, such that 


F(r,x) = - —Vlni?2(x) 


( 1 ) 


where r = and R^{x) is the product of the particle densities at the time of collision (r = 0). Note that 

this result is consistent with Ref. [s^, but the pre-factor is non-trivial (cf. @). The result ([T} has the added benefit 
that final particle momenta were found to be almost insensitive to the choice of the hydrodynamic starting time tq 
in heavy-ion collisions [s^. 

For the energy-density e(ro,x), initial conditions for each event are constructed as follows. Using Woods-Saxon 
distribution functions for the heavy ions such as An, Pb [s^ [s^ , the Hulthen wavefunction for the deuteron (cf. 0) 
and realistic calculations for the ^He wavefunction [^ . probability distributions of the nucleons within the nuclei 
of interest (cf. 0 ) are obtained. Using a Monte-Carlo Glauber 0 ], these probability distributions are mapped 
to positions of individual nucleons in the transverse (x, y) plane on an event-by-event basis implementing a hard¬ 
core repulsive potential of radius 0.4 fm between nucleons. The positions of nucleons undergoing at least one inelastic 
collisions are recorded (“participants”) and converted into a density function E? (x) by assuming that each participant 
contributes equally as a Gaussian with a width of re = 0.4 fm (to match the RMS radius of a single nucleon). This 
fixes the initial fluid velocities via Eq. o, and the energy density for hydrodynamics is assumed to be given as 

e('ro,x) = £'oi?^(x), (2) 


with Eq an overall constant (dependent on tq, collision energy and collision system) that is related to the total 
multiplicity of the event. Typically 100 event initial conditions are generated for each collision system. 

While the procedure to obtain the energy-density for hydrodynamics in Eq. Q is fairly standard in relativistic ion 
collision literature, the presence of the pre-equilibrium flow in Eq. © is relatively new. Eor this reason, also results 
without pre-equilibrium flow will be presented, and potential experimental ways to determine if pre-equilibrium flow is 
present in light-heavy ion collisions will be discussed. Other than the energy density and flow velocity field, relativistic 
hydrodynamic formulations also require the initial value for the shear and bulk stress tensors, which will be set to 
zero for simplicity. In heavy ion collisions, this assumption is harmless, as it has very little influence on final results 
cf. [H, However, given the short life-time of systems created in light-heavy ion collisions compared to heavy ion 
collisions, this assumption should be carefully revisited in follow-up studies. 


B. Thermal Stage: Viscous Hydrodynamics 

In modern definitions, hydrodynamics is understood to be an effective theory of energy conservation at long wave¬ 
length [H, H^l. Hydrodynamics is a good approximation of the bulk dynamics as long as higher order gradient 
corrections do not (strongly) change the leading order results. In standard nomenclature, viscous effects come into 
the energy-momentum tensor at first order in gradients (zeroth order would be ideal hydrodynamics). However, since 
all known consistent formulations for relativistic dissipative fluid dynamics are second order in gradients, this suggests 
a convenient handle on the effects of higher order gradients: the value of second-order transport coefficients, foremost 
the shear viscous relaxation time . 
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The view that is adopted here is the following: while first-order gradient effects in a system may be large (vis¬ 
cous effects sizable), a (viscous) hydrodynamic description of the system may still be quantitatively reliable as long 
as higher-order gradient corrections are small compared to first order gradients. This view is informed by direct 
simulations of strongly coupled quantum field theories out of equilibrium where it has been shown that viscous hydro¬ 
dynamics offers a reliable description even in regions where first-order (viscous) corrections to the ideal energy-stress 
tensor are approaching 100 percent [s^! 

To test if second-order corrections are small compared to first-order gradients, the value of Tjt (more precisely the 
ratio varied by 50 percent around the reference value C,, = 3, thereby generating a “systematic” error 

estimate of the applicability of viscous hydrodynamics. If the hydrodynamic gradient approximation was breaking 
down, one would expect second-order gradient terms to be as important as first-order gradient terms. Thus, final 
results for particle spectra should vary considerably when changing the strength of second order terms (via the value 
of Crj) by 50 percent. Conversely, if final results showed very little dependence on the value for Cr^ (as is the case 
for heavy ion collisions, cf. El), this could be considered evidence that a hydrodynamics was well applicable to such 
systems. In this sense, the reliability of hydrodynamics as a approximation to the system dynamics can be quantified 
and expressed as an error band generated by varying (7^, which is the strategy adopted in the following. 

Besides the second-order transport coefficients, the hydrodynamic evolution will also depend on choices for the 
(temperature-dependent) shear and bulk viscosity coefficients as well as the speed of sound (via the equation of 
state). Again, for simplicity constant values are used for the ratios of shear viscosity over entropy density rj/s and 
the bulk viscosity is set to identically zero (C = 0). For the equation of state, a parametrization of lattice QCD data, 

given in Ref. [i^, is employed. All of these choices should be revisited in a more detailed study. 

The hydrodynamic evolution is performed using the open-source code VH2-I-I (version 1.9) El|, adapted by smearing 
to prohibit code instabilities in the strong-gradient regions frequentl y en countered in event-by-event viscous hydrody¬ 
namics. The smearing (optimized from a version used before in Ref. [l0|) is performed by replacing low-energy density 
values by an average over nearest-neighbor cells, and it is only implemented at temperatures below the QCD phase 
transition (typically below 150 MeV). I have tested that final results are not sensitive to the details of the smearing 
implementation. The hydrodynamic evolution solves the equations of motion on a square lattice with area 20^ fm^ 
and 200^ grid points (lattice spacing O.I fm) and a time step of one hundredth of the lattice spacing. I have tested 

that the results are unchanged when simulating the same volume with a lattice spacing of 0.05 fm. 

During the hydrodynamic evolution, cells that cool below a certain switching temperature Tsw = 170 MeV are 
monitored. This condition defines a switching hyper-surface on which information about temperature, flow velocity, 
dissipative stress tensors as well as the normal vector of the hyper-surface are recorded for each cell. This information 
will provided the initial condition for the late-stage hadron cascade simulation. 


C. Late Stage: Hadron Interactions in a Hadron Cascade 

For the late-stage hadron interactions,the hadronic cascade code BSD E^ is used. Using the hyper-surface infor¬ 
mation to boost to the rest frame of each cell, the cascade is initialized with particles in the rest frame drawn from 
a Boltzmann distribution at a temperature Tsw with modifications of the momentum distribution to include defor¬ 
mations from viscous stress tensors (see E3 details). BSD includes hadron resonances in the particle data book 
up to masses of 2.2 GeV, which interact via simple s-wave scattering with a constant cross-section of 10 mb as well 
as scattering through resonances (modeled as a Breit-Wigner form). Once the resonances have stopped interacting, 
one can obtain final charged hadron multiplicities mean charged particle momentum (px) and flow coefficients 

Vnipx) for n>l from summing over individual particles with momenta p. Specifically, 

E ch. particles ^oo dN^ 

_ _ in PT bin dNc\i _ dNc\i /-nr^\ _ Jq dYdpj^ 

27TpTdYdpT 27 TPt^T^Y 

\Vn\{PT) = \/Sn(PT)^ + C„(pt)^ , 

where At = 80 MeV, Ay = 2 are the width of bins for particle px and rapidity V, respectively. Note that since 
the cascade is applied to a boost-invariance case, the large Ay value is of no significance. In practice, a sum over 
both particles and anti-particles and division of the spectra by two is performed, in order to increase statistics. For 
every hydrodynamic evolution event, 100,000 BSD events are run to increase statistics. In doing so, the sums in the 
definition of Vn above are extended over all BSD events, thereby explicitly ignoring fluctuations arising from hadronic 
decays. After thus obtaining results Vn{px) for each hydrodynamic event, an event average to obtain 


, dY ~ JO dYdpT ’ 


SniPr) \ ^ 

Cnipr) J 


E 


ch. particles 
p-p bin 


sin(n(/))) 

cos(n^) 


e: 


ch. particles 
prp bin 


-JY- 


' = arctan 


(£) 


( 3 ) 
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System yi [GeV] o-jvAr[mb] (A''part) rijs preflow? ^^(th) _ ^^(exp) _ (pr) [GeV] (th) (pr) [GeV] (exp) 


p-fAu 

7.7 

32 

7.4 

0.08 

yes 

4.5 


0.550(6) 



p-fAu 

7.7 

32 

7.4 

0.08 

no 

4.6 


0.497(8) 



p-fAu 

62.4 

36 

8.6 

0.08 

yes 

8 


0.588(5) 



p-fAu 

62.4 

36 

8.6 

0.08 

no 

8.4 


0.528(6) 



p-fAl 

200 

42 

4.7 

0.08 

yes 

5.4 


0.575(6) 



p-fAl 

200 

42 

4.7 

0.08 

no 

5.5 


0.515(7) 



p-fAu 

200 

42 

9.6 

0.08 

yes 

10 


0.604(5) 



p-fAu 

200 

42 

9.6 

0.08 

no 

10 


0.537(5) 



d-l-Au 

62.4 

36 

14.9 

0.08 

yes 

15 


0.570 



d-l-Au 

200 

42 

17.5 

0.08 

yes 

20 

20.3 ± 1.7 (0-5%)[llJ 

0.576(1) 

0.554 

(0-20%) [4^ 

d-l-Au 

200 

42 

17.5 

0.08 

no 

20 

20.3 ± 1.7 (0-5%) [UJ 

0.523(5) 

0.554 

(0-20%) [4^ 

®He-tAu 

62.4 

36 

21 

0.08 

yes 

21 


0.557(1) 



^He-tAu 

62.4 

36 

21 

0.08 

no 

21 


0.509(4) 



^He-tAu 

200 

42 

24 

0.08 

yes 

27 ± 1 


0.567(2) 



®He-tAu 

200 

42 

24 

0.08 

no 

27 ± 1 


0.520(3) 



p-fPb 

5020 

70 

15.2 

0.16 

yes 

39 ± 1 

35 ±0.5 (0-20%) [1] 

0.716(6) 



p-fPb 

5020 

70 

15.2 

0.08 

no 

38 ± 1 

35 ±0.5 (0-20%) [IJ 

0.623(3) 



®He-tPb 

5020 

70 

32.4 

0.16 

yes 

74 


0.676 



TABLE I. Comparison of superSONIC 

runs (“th”) with available experimental data (“ 

exp”). For all model runs w = 0.4 fm. 

To = 0.25 fm/c and Tsw 

= 0.17 GeV as 

well as 

“central collisions” with impact parameter b < 2 fm/c 

are chosen. The effect 


of changing C,, = 2 — 3 is contained in the reported theoretical error estimates. Mean particle pr is for unidentified charged 
hadrons. For the experimental data both the centrality and the reference are reported. 


the event-by-event mean and event-by-event fluctuation is performed, the latter of which is recast into a statistical 
error bar on the mean. The results from this procedure are reported on in the following. 

The combined simulation package of the early-stage, thermal stage and late stage evolution thus described above 
will be referred to as “superSONIC” in the remainder of this article. 


III. RESULTS 

A summary of systems that were simulated is given in Tab. ID The first column in this table gives the system 
configuration, the second the collision energy, the third the inelastic cross-section (from Ref. 0) used in the Monte- 
Carlo Glauber. All Monte-Carlo Glauber events were generated for “central” collisions by imposing an impact 
parameter b < 2 fm, loosely corresponding to the 0-5% most central collisions. The fourth column in Tab. |T] refers 
to the mean number of participants obtained by averaging over 100 random Monte-Carlo Glauber events. The fifth 
column specifies the constant value of p/s chosen in the hydrodynamic simulations, which were all started at an 
initial time tq = 0.25 fm/c with or without pre-equilibrium flow according to the sixth column of Tab. HI For all 
runs, the energy scale factor Eq that was chosen in order to match measured or expected values ([1^) for the charged 
particle multiplicity, reported in column seven and eighth of Tab. HI (Note that the calculated is converted 

to experimentally measured pseudo-rapidity distribution by dividing by 1.1). For experimentally measured 

quantities, also the centrality class for the measurement as well as the corresponding reference is reported. Finally, 
the last two columns in Tab. HI give the mean charged particle transverse momentum in the superSONIC simulation 
compared to experimental data where available. 

As one can see from Tab. HI simulated particle multiplicities range from ~ 74 down to ~ 4.5. This implies 
that systems with very few particles are being simulated and one generally expects hydrodynamics to be less applicable 
to these fewer-body systems. Also, note that there is a clear change in the mean charged particle momentum for 
systems with compared to without pre-equilibrium flow. This is not too surprising given that systems created in light- 
heavy ion collisions live comparatively shorter than those created in heavy ion collisions, thus making light-heavy ion 
collision systems more sensitive to pre-equilibrium conditions. While other parameters (viscosity, choice of switching 
temperature T^w) also affect the particle mean momentum, the strong difference between results with and without 
pre-equilibrium flow could serve as important discriminatory tool that is easy to implement experimentally. 

In figured! results for the flow coefficients Vn, n=2,3,4 minus are shown for p-|-Pband ^He-|-Pbcollisions at 
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•y/s = 5.02 TeV per nucleon pair (LHC energies). Results are reported as a difference with respect to calculated 
(the highest harmonic calculated) because finite statistics from the 100,000 B3D runs start to pollute the results at 
high pt and the true v^(j)t) can reasonably be expected to be consistent with zero. Hence the calculated v^{pt) 
is a good measure of the statistical error, and can be used to subtract the statistical fluctuation for the other flow 
harmonics. In principle, this procedure could be made unnecessary by rerunning all the simulations with at least 10® 
B3D runs per hydro event, which is left for future work. 

For comparison, results with and without pre-equilibrium flow are shown. Boxes indicate uncertainty arising from 
both statistical fluctuations as well as systematic errors, the latter of which are quantified by performing simulations 
at different values of C,; = 2 — 3. From these plots, the first finding is that the hydrodynamic uncertainty range 
for light-heavy ion collisions at ^/s = 5.02 TeV is rather small compared to the overall magnitude of the flow for 
V 2 ,V 3 -,V 4 , thus providing solid theoretical support to the notion of true hydrodynamic behavior in light-heavy ion 
collisions. Furthermore, one finds that simulated flow coefficients are in overall agreement with experimental data, 
where available, when suitably adjusting the constant simulated shear viscosity over entropy density. Note that, 
somewhat surprisingly, it is hard to distinguish the cases of 77 /s = 0.16 with pre-equilibrium flow and 77 /s = 0.08 
without pre-equilibrium flow using experimental data for 772 , 773 , 774 - For this reason alone, it would be highly desirable 
to have an experimental handle on the presence of pre-equilibrium flow. Finally, comparing simulations of ®He-|-Pb and 
p-fPb collisions at LHC energies one finds that one could expect V 2 to be 50 percent higher in ®He-|-Pb collisions than 
in p-fPbcollisions, while the overall magnitude for 773,774 would be approximately the same. This presumably points 
to the fact that in ®He-|-Pb collision there is a sizable geometric component (other than event-by-event fluctuations) 
that drives V 2 - 

In figure [3l results for the flow coefficients 77 „, n=2,3,4 minus 775 are shown for p-fAl, p-|-Au, d-l-Auand 
®He-|-Aucollisions at -/i = 200 GeV per nucleon pair (highest RHIC energies). As before, results with and without 
pre-equilibrium flow are shown and boxes indicate theoretical uncertainty. Again, one finds that simulated flow 
coefficients are in overall agreement with experimental data, where available. Interestingly, the simulated results, 
both with and without pre-equilibrium flow predict a 773 component in d-|-Au collisions that is almost as large as 
in ®He-|-Au. While 773 in central ®He-|-Auhas been measured by the PHENIX experiment (@), a non-vanishing 
773 in d-|-Au central collisions has not been seen in any experiment yet. However, given the unambiguous presence 
of a sizable 773 component in the superSONIC simulation package, as well as in other theory (hydrodynamic and 
non-hydrodynamic) simulations [l3l Ifilj , measuring or putting an upper bound on 773 in d-f Au collisions at ^ = 200 
GeV could serve as a very useful experimental verification of theory ingredients to light-heavy ion collisions. 

Similar to the finding for ^ = 5.02 TeV, for results at y/s = 200 shown in Fig. |3]one finds noticeably smaller V 2 
component in p-fAland p-|-Aucollisions compared to d-|-Auand ®He-|-Au. The comparison between different systems 
for 773 is highly dependent on the presence or absence of pre-equilibrium flow. For instance, with pre-equilibrium 
flow, one finds that 773 in p-|-Au collisions at = 200 GeV is almost as large as for ®He-|-Au, while without the 
presence of pre-equilibrium flow 773 in p-|-Auis about half as large as in ®He-|-Au. This suggests that 773 in light-heavy 
ion collisions at y/s = 200 GeV could be a good experimental handle on pre-equilibrium QGD dynamics. Unlike the 
situation found for ^/s = 5.02 TeV, simulation results predict 774 to be consistent with zero except for p-|-Au ,d-|-Auand 
®He-|-Au where some small, non-vanishing 774 is found. Since 774 is found to be so small at ^ < 200 GeV, it is likely 
to be hard to measure, and hence it will be discounted as a probe for pre-equilibrium dynamics in the following. 

What is striking about the results shown in Fig. [3] is the magnitude of V 2 (and even 773 in the case of pre-equilibrium 
flow) predicted in p-|-Al collisions at y/s = 200 GeV. The system created in these collisions consists of an average 
multiplicity of only ^ 5.4, yet superSONIG results exhibit a clear flow response much larger than the estimated 
uncertainty band for the applicability of hydrodynamics. According to the criterion defined above, I predict that 
systems created in central p-|-Al collisions at y/s = 200 GeV can be quantitatively described using viscous hydrody¬ 
namics, and that a V 2 flow component in these systems can be expected to be similar to that for p-|-Au collisions at 
Vs = 200 GeV. 

In Fig. m results for the flow coefficients 77„, n=2,3,4 minus 775 are shown for p-|-Au, d-|-Auand ®He-|-Aucollisions at 
= 62.4 GeV per nucleon pair. Gomparing results in Fig. Uto results at ^/s = 200 GeV in Fig. |31 it is hard to find 
clear differences in any of the flow observables shown. All 77„ at ^/s = 62.4 GeV in all systems are broadly consistent 
with results at ^/s = 200 GeV, only slightly lower. Nevertheless, one should point out that this in particular implies 
a sizable V 2 component in p-|-Auat ^/s = 62.4, and hence hydrodynamic behavior at these collision energies. 

In order to answer the question of a breakdown of hydrodynamic applicability, results for flow coefficients in 
proton-nucleus collisions (p-|-Auand p-|-Pb) at energies ranging from ^/s = 7.7 GeV to ^/s = 5.02 TeV are compared 
in figure [5] Most remarkable, the simulations predict a sizable 772 ( 2.5 px) ~ 5 — 8 percent even for p-|-Aucollisions 
at y/s = 7.7 GeV, with a predicted multiplicity of only ~ 4.5. Again, according to the criterion of applicability 
of hydrodynamics from above, hydrodynamics still applies for systems created in p-|-Au collisions at y/s = 7.7 GeV. 
Failure of finding a break-down of hydrodynamics in p-|-Au collisions at y/s = 7.7 GeV begs the question of which 
value of y/s, if any, one would have to study in order to truly see hydrodynamic no longer apply. Unfortunately, at 
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FIG. 1. Temperature evolution in systems created in light-heavy ion collisions. Shown are temperature values encountered in 
event-by-event superSONIC simulations in ®He-|-Au ,®He-|-Pb, p-|-Auand p-|-Pb collisions at various energies (bars) and event- 
by-event mean starting (maximal) temperatures < Tstart > at r = 0.25 fm/c. Note that p-|-Au collisions at RHIC energies 
probe a smaller range of temperatures, thereby being insensitive to transport properties for T > 0.33 GeV. 


collision energies below y/s = 7.7 GeV, the zero chemical potential lattice equation of state employed in superSONIC 
is clearly no longer applicable, so studying proton-nucleus collisions at energies below y/s = 7.7 GeV is not feasible 
with the current approach. 

While a complete break-down of hydrodynamics is not observed in superSONIC simulations, one does observe a 
gradual break-down of hydrodynamics. That is, starting with p-|-Pb collisions at y/s = 5.02 TeV and lowering the 
collision energy down to y/s = 7.7 GeV for p-|-Au collisions, first Vi and then first drop and then collapse to values 
consistent with zero for all momenta considered. In the same fashion, one does also observe a clear drop in V 2 as a 
function of lowering y/s, even though at y/s = 7.7 GeV V2 results have not (yet) collapsed to zero. Hydrodynamics is 
breaking down, but is has not broken down completely for p-|-Au collisions at y/s = 7.7 GeV. My interpretation of this 
finding is that the question of applicability of hydrodynamics or collectivity in small systems does not have a yes/no 
answer, but rather should be thought of as a gradual process similar to the confinement/deconfinement cross-over 
transition in QCD, where the value of the critical temperature is also dependent on the observable one considers. 

In Fig. ini the flow response V 2 , vs is compared for different collision systems in an attempt to quantify the presence 
of pre-equilibrium flow. Similar to the results shown in figures above, one finds that V 2 is fairly insensitive to pre¬ 
equilibrium flow, regardless of the collision energy. However, turns out to be a good indicator for the presence 
of pre-equilibrium flow in different collision systems at y/s < 200 GeV. Specifically, while V 3 in ^He-|-Au collisions 
has been measured at y/s = 200 GeV, a measurement of V 3 in d-|-Au collisions at that same energy could serve 
as an experimental handle on the presence of pre-equilibrium flow. Also, note that lower collision energies (such 
as y/s = 62.4 GeV) would be even better suited for an experimental test of pre-equilibrium flow since the systems 
created in these collisions live comparatively shorter, and are thus more sensitive to out-of-equilibrium QGD transport 
dynamics. Gonversely, it would be harder to probe these pre-equilibrium transport effects at y/s = 5.04 TeV, both 
because ^He-I-Pb collisions at the LHC are not currently planned, and because ^He-I-Pbto p-l-Pbua ratios are only 
sensitive to pre-equilibrium flow at px ^ 0.5 GeV. 

Another point worth noting is that the more limited temperature range encountered in light-heavy ion collisions as 
compared to heavy-ion collisions, especially at RHIC energies. Fig. [1] shows the maximal temperature encountered in 
various collision systems from y/s = 7.7 GeV to y/s = 5.02 TeV. Also shown in Fig. [1] is the typical (event-by-event 
mean) starting temperature, which is a measure of the maximal temperature the average system starts out with 
in these collisions (dots). One finds that hydrodynamic evolution in p-|-Pbcollisions probes and thus averages over 
transport properties (such as rj/s and (/s) over a large temperature range T < 0.53 GeV (< T >= 0.4 GeV, while by 
contrast p-|-Aucollisions at y/s = 7.7 GeV only probe T < 0.33 GeV (< T >= 0.25 GeV). Hence p-|-Aucollisions at 









^/s = 7.7 GeV are mostly sensitive to QCD transport properties at T < 0.25 GeV, which could be a key experimental 
handle on probing the temperature dependence on e.g. ry/s. 

Finally, results on HBT Radii for selected collision systems are reported in Tab. H] in Figures [7l8] and El 


IV. SUMMARY AND CONCLUSIONS 

In this work, central light on heavy ion collisions from ^ = 7.7 GeV to ^/s = 5.02 TeV were simulated using 
superSONIC, a model that combines pre-equilibrium flow, relativistic viscous hydrodynamics and a hadron cascade 
afterburner. By varying the strength of the second order transport coefficients, one could quantify to which extent 
viscous hydrodynamics offers a reliable description of these systems . It was found that even in p-|-Au collisions at 
y/s = 7.7 GeV, a sizable collective flow component V 2 much larger than the systematic hydrodynamic uncertainty 
is present. Thus, viscous hydrodynamics is still applicable to describe V 2 in the systems created in these small, 
low-energy, few-body collisions. However, there is evidence that hydrodynamics is breaking down gradually as y/s is 
lowered from 5.02 TeV to y/s = 7.7 GeV. Specifically, first U 4 , then V 3 and finally V 2 start to decrease and eventually 
V 3 and Vi become consistent with zero as the collision energy is lowered. The question of whether hydrodynamics 
applies to describe light-heavy ion collisions therefore cannot be answered by a simple yes or no, but depends on 
the quantity in question. In future work, it would be interesting to perform simulations at collision energies below 
yi=7.7 in order to potentially observe also V 2 collapse to zero in p-|-Au collisions. 

For d-l-Au collisions at y/s = 200 GeV superSONIC simulations predicted V 3 results that were of the same order 
of magnitude as those measured in ^He-|-Au. Since similar results are observed in other theoretical models for 
d-|-Au collisions at y/s = 200 GeV, I predicted that V 3 should be observable in experiment. Furthermore, the sensitivity 
of flow results to the presence of pre-equilibrium flow was studied. It was found that low-energy collisions were most 
sensitive to the presence of pre-equilibrium flow and pointed out that the V 3 ratio of ^He-|-Auto d-|-Au could provide 
an experimental handle on pre-equilibrium QCD dynamics. 

Moreover, simulation results clearly show that light-heavy ion collision systems probe a temperature window focused 
around the QCD phase transition temperature. Thus, a combination of simulation results and experimental data for 
these systems would offer a promising handle on the temperature dependence on QCD transport properties. 

Many aspects of this study are amenable to improvements in future work, but despite the current set of limitations 
and approximations, the findings of this work could hopefully demonstrate the usefulness of studying light-heavy ion 
collisions to learn about QCD. 
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Appendix A: Plot Collection 

For better readability of the text, many of the plots for this article are collected in this appendix. 

LHC, = 5.02 TeV 


superSONIC: p+Pb @ 5.02 TeV, t|/s=0.1 6 superSONIC: pPb @ 5.02 TeV, ti/s=0.08 






FIG. 2. Flow harmonics Vnipr) — vsipr) for unidentified charged particles (unid) for n = 2, 3,4 from superSONIC, with and 
without pre-equilibrium flow. Boxes indicate combined statistic and estimated systematic error for hydrodynamics (latter from 
varying Crj = 2 — 3). For reference, experimental data is shown where available am- 
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RHIC, = 200 GeV 
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superSONIC: p+AI (§) 200 GeV, ti/s=0.08 





0 0.5 1 1.5 2 2.5 

pt [GeV] 


superSONIC: ^He+Au @ 200 GeV, ti/s=0.08 




0 0.5 1 1.5 2 2.5 

pt [GeV] 


superSONIC: p+Au @ 200 GeV, ti/s=0.08 





pt [GeV] 


FIG. 3. Flow harmonics Vn(pT) — V 5 (pt) for identified pions (tt) and unidentified charged particles (unid) for n = 2,3,4 from 
superSONIC, with and without pre-equilibrium flow. Boxes indicate combined statistic and estimated systematic error for 
hydrodynamics (latter from varying C-q = 2 — S). For reference, experimental data is shown where available [3, . 
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RHIC, = 62.4 GeV 



0 0.5 1 1.5 2 2.5 

pt [GeV] 


superSONIC: d+Au @ 62.4 GeV, ti/s=0.08 




superSONIC; p+Au @ 200 GeV, ti/s=0.08 





FIG. 4. Flow harmonics v„{pt) — V5{pt) for identified pions (tt) and unidentified charged particles (unid) for n = 2,3,4 from 
superSONIC, with and without pre-equilibrium flow. Boxes indicate combined statistic and estimated systematic error for 
hydrodynamics (latter from varying Cn = 2 — 3). 



















































































































































































































































Proton-Nucleus collisions, yi=7.7-5000 GeV 
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superSONIC: v2 for p+A @ 7.7 to 5020 GeV 


superSONIC: v2 for p+A @ 7.7 to 5020 GeV 
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superSONIC: v3 for p+A @ 7.7 to 5020 GeV 


superSONIC: v3 for p+A @ 7.7 to 5020 GeV 
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superSONIC: v4 for p+A @ 7.7 to 5020 GeV 


superSONIC: v4 for p+A @ 7.7 to 5020 GeV 
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FIG. 5. Flow harmonics Vnipr) — v^(pt) for unidentified charged particles (unid) for n = 2, 3,4 from superSONIC in proton- 
nucleus collisions, with and without pre-equilibrium flow. Boxes indicate combined statistic and estimated systematic error for 
hydrodynamics (latter from varying G,, = 2 — 3). Y-axis is same scale on all plots. The p+Pb system, which has the highest 
multiplicity (see. Tab.[T|, shows non-vanishing flow components up to n = 4, but as the multiplicity is decreased, first U 4 , then 
Vi and eventually also V2 start to decrease and (in the case of 114) eventually become consistent with zero. 
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superSONIC: p,cl,^He+Au @ 62.4 GeV 


superSONIC: p,cl,^He+Au @ 62.4 GeV 
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superSONIC: p,He®+Pb @ 5.02 TeV 
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FIG. 6 . Ratios of V 2 ,V 3 between ®He+Au,p+Auand d+Au(and ®He+Pband p+Pb), with and without pre-equilibrium 
flow (“preflow”). Boxes indicate combined statistic and estimated systematic error for hydrodynamics (latter from changing 
C 77 = 2 — 3). For < 200 GeV, the 1/3 ratios between ^He-I-Auand p+Auto d+Auare rather sensitive to the presence of 
preflow. For ^/s = 5.02 TeV, this difference is much smaller since the evolution is less sensitive to the presence of pre-equilibrium 
flow. 





















































































































































































































































































































































































































































































LHC, = 5.02 TeV 
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superSONIC: p+Pb (2) 5.02 TeV, ti/s=0.16 
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FIG. 7. HBT Radii for pions for p+Pb and ^He+Pb collisions at ^/s — 5.02 TeV. Varying Ctj = 2 — 3 in hydrodynamics results 
in changes smaller than the symbol size shown. 
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superSONIC: p+Au @ 200 GeV, ti/s=0.08 
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FIG. 8. HBT Radii for pions for p+Au, d+Auand ^He+Au collisions at y/s = 200 GeV. Waring Gt? = 2 — 3 in hydrodynamics 
results in changes smaller than the symbol size shown. 
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superSONIC: d+Au @ 62.4 GeV, t]/s=0.08 
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FIG. 9. HBT Radii for pions for d+Au and ®He+An collisions at y/s = 62.4 GeV. Waring C,, = 2 — 3 in hydrodynamics results 
in changes smaller than the symbol size shown. 










